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Abstract: We study the hydrodynamic properties of strongly coupled SU(N) 
Yang- Mills theory of the Dl-brane at finite temperature in the framework of 
gauge/gravity duality. The only non-trivial viscous transport coefficient in 1 + 1 
dimensions is the bulk viscosity. We evaluate the bulk viscosity by isolating the 
quasi-normal mode corresponding to the sound channel for the gravitational back- 
ground of the Dl-brane. We find that the ratio of the bulk viscosity to the entropy 
density to be l/Air. This ratio continues to be 1/Att also in the regime when the Dl- 
brane Yang-Mills theory is dual to the gravitational background of the fundamental 
string. Our analysis shows that this ratio is equal to 1/4-7T for a class of gravita- 
tional backgrounds dual to field theories in 1 + 1 dimensions obtained by considering 
Dl-branes at cones over Sasaki-Einstein 7-manifolds. 
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1. Introduction 

Recent studies of gauge theories using the framework of the AdS/CFT correspon- 
dence has revealed that gauge theories at large 't Hooft coupling and at long-distances 
and low- frequencies can be described by fluid mechanics |0 3 @) @> B S @> 01 1 - ^ 
hydrodynamic description implies that the correlation functions of components of 
stress-energy tensor or conserved currents are fixed once a few transport coefficients 
are known. These transport coefficients have been evaluated for several examples 
of gauge theories at strong coupling using the AdS/CFT correspondence. Recently, 

1 For a review and a complete list of references please see [p| 
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transport coefficients for non-linear hydrodynamics have been obtained || |10[. These 
studies indicate that for field theories which admit a gravity dual, the ratio of shear 
viscosity 77 to the entropy density s at strong t'Hooft coupling is universal |5|, [TT] and 
is given by 

- — j-. (i.i) 

Most studies of this ratio have been focused on asymptotic anti-de Sitter backgrounds 
in various dimensions. The gauge/gravity correspondence also applies to stacks of 
Dp-branes for arbitrary p |T2|, |13[ 2 . It was shown in [R, [L5| that the ratio of shear 



viscosity to entropy density for theories on Dp-branes for p > 2 also is 1/4tt 3 . An 
investigation of transport coefficients for the theory on the Dl-brane is missing in 
the literature. 

In this paper, we begin a study of the hydrodynamic behaviour of field theories 
in 1 + 1 dimensions which admit a gravity dual. What makes one spatial dimension 
special is the absence of shear. In fact, for the conformal field theories in 1 + 1 di- 
mensions, there are no transport coefficients. This is because there are no non-trivial 
components of the symmetric traceless stress-energy tensor in these dimensions. The 
stress tensor is that of a perfect fluid. Thus, to study non-trivial transport proper- 
ties of field theories in 1 + 1 dimensions, it is necessary to study the non-conformal 
case. The only viscous transport coefficient for non-conformal field theories in 1 + 1 
dimensions is the bulk viscosity. The simplest example of such a non-conformal field 
theory is the theory on the Dl-brane, the 1 + 1 dimensional SU(N) gauge theory 
with 16 supersymmetries. It can be obtained as a dimensional reduction of M = 4 
SYM from 3 + 1 dimensions. We consider this theory at finite temperature. It admits 
dual gravity descriptions in two regimes 

(i) V\N~ 2/3 « T « and (1.2) 

(ii) VXN- 1 « T « V\N~ 2/S . 

Here, the 't Hooft coupling is denoted by A = gy M N and T is the temperature. In 
regime (i), the gravity dual is the background of that of the non-extremal Dl-brane, 
while in regime (ii), the dual is that of non-extremal fundamental string. 

One reason why transport coefficients were not studied for backgrounds corre- 
sponding to the Dl-brane is that the gauge invariant fluctuations of supergravity 
fields found for the case of Dp brane with p > 2 cannot be extended to p — 1. In 
this paper, we isolate the gauge invariant fluctuation of the graviton and the dilaton 
which corresponds to the sound channel. From its equation of motion, we find the 
following dispersion relation for its quasinormal mode 

oo = ^-^-q 2 - (i-3) 



2 For a review and recent developments on holography for the non-conformal case, please sec [14 
3 Other non-conformal systems were studied in |l6|, [l?], [l8| 



-2- 



This dispersion relation for the quasinormal mode remains the same for both the 
Dl-brane background as well as for the Fl-string background. We show that the 
retarded two point functions of components of stress tensor has a pole at the value of 
uj corresponding to the sound mode. From the universal properties of hydrodynamics 
in 1 + 1 dimension, we find that the pole is given by 

lu = v s q — i , - — -q 2 . (1.4) 

Here, v s is the speed of sound in the medium, e the energy density and P the pressure. 
By comparing the dispersion relation ( |1.3| ) and the above equation, we can read out 
the following properties of the Dl-brane gauge theory in both the regimes given in 

(ED- 

Vs = ^ \ = h' (L5) 

where £ is the bulk viscosity and s is the entropy density. As a cross check of our 
calculations we use the Kubo's formula to evaluate the bluk viscosity and show that 
the ratio £/s is given by 1/4tt. 

It is curious that the ratio of the bulk viscosity to entropy density has the same 
value as that of ratio rj/s for strongly coupled field theories which admit gravity 
duals in higher dimensions. Also the fact that the ratio £/s remains 1/Atc for both 
the Dl-brane and the Fl-string suggests that this ratio might be universal for a class 
of gravitational backgrounds. We show that the this ratio continues to be 1/471 for 
the class of 1 + 1 dimensional non-conformal field theories on Dl-branes at cones over 
Sasaki-Einstein 7-manifolds. 

The organization of the paper is as follows: In the next section, we briefly review 
the gauge/gravity correspondence and thermodynamics for the case of the Dl-brane. 
In section 3, we discuss the implications of hydrodynamics on the thermal Green's 
functions in 1 + 1 dimensions. In section 4, we discuss the details of how to isolate 
the gauge invariant fluctuation of the dilaton and the graviton corresponding to the 
sound channel. We then derive the dispersion relation of the quasi normal mode for 
the sound channel. We also evaluate the two point functions of the stress tensor 
components from gravity and show that the Lorentz structure and the structure of 
its pole agrees with that predicted from general considerations of hydrodynamics in 
1 + 1 dimensions. In section 5, we show that the ratio of £/s is l/4ir for both the Dl- 
brane background and that of the fundamental string. This is done in two ways: we 
match the dissipative part of the quasi-normal mode for the sound channel in gravity 
to that expected from general hydrodynamic considerations and more directly by 
applying the Kubo formula for bulk viscosity in terms of stress tensor correlators. 
Both methods yield the ratio of bulk viscosity to entropy density to be 1/Att. We 
then show that this ratio continues to be 1/47T for the class of field theories in 1 + 1 
dimensions dual to Dl-branes at cones over Sasaki-Einstein manifolds. Appendix A. 
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shows that the constraints we impose on the graviton and the dilaton perturbations of 
the Dl-brane background are consistent with their equations of motion. Appendix B. 
contains proofs of the two identitites which are used in our derivation of the equation 
for the sound mode. 



2. Gauge/gravity duality for the Dl-brane 

In this section, we briefly review the statement of gauge/gravity duality for the case 
of N Dl-branes. This enables us to set the notations and conventions that we will use 
and also to state the bounds on temperatures for which the bulk viscosity evaluated 
holographically from the gravity background can be trusted. 

In it was argued that SU(N) Yang-Mills with 16 supercharges in 1 + 1 
dimensions at large N is dual to the near horizon supergravity solution of Dl-branes. 
The near horizon supergravity solution in the Einstein frame is given by 



ds 2 

P 4>(r) 



*F: 



RR 



H^{r){-dt 2 + dx\) + m{r){dr 2 + r 2 d{l 2 ) 



(2.1) 



where 



H(r) 



and V = g^ M 2 6 7r 3 Na 



./i 



(2.2) 



and dfl 2 refers to the metric on the unit 7-sphere and u>s 7 its volume form. Note 
that, we have dualised here the Ramond-Ramond charge of the Dl-brane. The 
gravity description is valid in the energy domain 



9y M N* « U « qym^/N- 



(2.3) 



Here U sets the energy scale. Note that the Yang-Mills coupling in 1 + 1 dimensions 
has the units of energy. For completeness, we mention that the background in ( [2.1Q is 
the solution of type IIB supergravity equations of motion in 10 dimensions obtained 
from the following action 



S 



1 



IIB 



167TG10 



-g 



R {9) ~ 



2-3 



■ e t>{F^ R) ) 2 



(2.4) 



At both ends of the domain in Q2.3p , the curvatures of the supergravity solution in 
( |2.1|) grow and the solution breaks down. In the UV i.e. U » gyMvN, one can 
trust the perturbative description of Yang-Mills theory at any value of N. However 
in the domain 



gvu « U « g Y MN'( 



(2.5) 
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the dual description is given by the near horizon geometry of the fundamental string 
solution. This background is obtained by performing a S-duality transformation on 
the D-brane background in (|2.1| ) 

ds 2 10 = H--±(r)(-dt 2 + dx\) + H*(r)(dr 2 + r 2 dn 2 ), (2.6) 
e Hr) = H(r)-?, 
*F 3 NS = 6L 6 u St . 

Note that the only changes in the background compared to that of the Dl-brane is 
cf) — > —0 and the Ramond-Ramond flux replaced by the Neveu-Schwarz flux on the 
7-sphere. The above background is a solution to the equations of motion from the 
following action 



1 ' ,io. 



Sub = , a n / d x^g 

107rGio 



R{g)-\d M <pd M <P-^e~\Ft 



(2.7) 

Unlike the case of the AdS / CFT duality, the supergravity solution in (|2.1|) and p] 



is not asymptotically AdS3 but only conformal to AdS3 fEzf . Finally, deep in the 
IR i.e. U « Qym i the valid description is given in terms of the conformal field 
theory on the orbifold (R 8 ) N /Sn for any N. Thus both in the UV and in the IR, 
the 1 + 1 dimensional super Yang-Mills flows to a conformal field theory. It is only 
in the domain qym « U « gyuvN and in the limit of large N, one has a dual 
description in terms of a supergravity solution. 

To study hydrodynamics of the Dl-brane theory, we need to consider the theory 
at finite temperature. The dual description is given in terms of the near horizon 
geometry of the non-extremal Dl-brane solution which is given by 

dsj = H--*(r)(-f(r)dt 2 + dx{) + H*(r) (^ + r 2 d{l 2 ^ , (2.8) 

e 0(O = H(r)K 
F 7 RR = 6L 6 u S7 , 

where 

m = i - (7) ■ ( 2 -9) 

The temperature and the entropy density of the Dl-brane theory is related to the 
non-extremal parameter tq by 

T- 3r ° 2 a- 1 T°V- ^ r 4 L 3 (2 10) 

We can convert the domain of the validity of the Dl-brane solution to ranges in 
temperature by identifying Uq = r /a' with the U in Q2.3D, using the defintion of L 
in ( p.2[) and the following relations. 

9ym = 2^7' G10 = 2 vr g s a . (2.11) 
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We obtain 

V\N^ « T « (2.12) 

Here, we have defined the t'Hooft coupling A = (/y M iV. It is now clear that for large 
N, this is a sufficiently large domain. Now for T << -\/AiV _2 / 3 , the holographic dual 
of the Yang-Mills theory is given by the non-extremal fundamental string solution. 

dsi = H-l(r)(-f(r)dt 2 + dxl) + tfi(r) + rW^j , (2.13) 

e 4>(r) = H(r)-^, 
F 7 NS = 6L 6 ^ 7 . 

Again writing the domain in (|2.5| ) in terms of temperature and the 't Hooft coupling, 
we obtain that the above solution can be trusted in the following temperature range 

VXN- 1 « T « VXN-I . (2.14) 

To conclude, for very high temperatures T » v^A and for very low temperatures 
T « the Yang-Mills theory flows to a free conformal field theory. 

In this paper we evaluate the bulk viscosity £ using the holographic description 
for temperatures in the regime « T « v^A- From the above discussion, 

we see that for T » \/~X and T « V% the theory flows to a conformal field theory. 
Therefore, we expect the bulk viscosity to vanish in these domains. However, in the 
regimes (|2.12| ) and (|2.14|) , we will see that we obtain a non-trivial value of the bulk 
viscosity. 



3. Hydrodynamics and the sound mode in 1 + 1 dimensions 

In this section, we will discuss generalities of relativistic hydrodynamics in 1 + 1 
dimensions. As mentioned in the introduction, hydrodynamics in 1 + 1 dimensions is 
special due to the absence of shear. We discuss here the constraints of conservation 
laws and hydrodynamics on the structure of the thermal Green's function of the stress 
tensor in 1 + 1 dimensions. We show that this implies that the only hydrodynamic 
mode is longitudinal and we will determine its dispersion relation. 

3.1 Lorentz structure of the correlators 

We show using translational invariance and conservation of the stress tensor that its 
Green's function is entirely characterized by a single function in 1 + 1 dimensions. 
We define the retarded Green's function of the stress energy tensor to be 

(WO* -V) = ~ i9 ( x ° ~ y°)([T^),T aP (y)]). (3.1) 

Making use of translation invariance of the state, we can define the Fourier trans- 
form of the above correlator denoted as G ^^(k) . It is symmetric by definition in 
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indices (/x, v) and (a, (3). Further more, we have the following symmetry due to CPT 
invariance. 

(3.2) 

The conservation of the stress-energy tensor leads to the following Ward identity 

WG^ aP {k) = 0. (3.3) 
This suggests a useful tensor which forms a basis to write down the correlator is 

P*, = ^ - ^r- (3-4) 

Note that k^P^ u = 0. If the states involved in the expectation value in ( |3.1|) are 
Lorentz invariant, a convenient decomposition is to split the components into a part 
which contains the trace Tj^rj^G^ap and the traceless. This is given by 



where 



G^ a p{k) — Pfj,uPapGB(k ) + H^ a pGs{k ), (3.5) 

Note that rj^H^^p = and the two tensors in ( |3.5| ) are orthogonal; P^PapH^ *, = 
0. At this stage, it seems neccessary that one needs 2 functions to characterize the 
2-point function of the stress energy tensor namely, Gb and Gs- Now substituting 
explicitly the value of k^ = (—u, q), we find that 

P - g 2 P - ^ P - ^ (* 7 ) 

uj z — q z oo z — q z oo z — q z 

Using the above components of the tensor P^, it is easy to see that all components 
of H^ iu a p vanish. Therefore, the two point function of the stress tensor in a 1 + 1 
dimensional theory is entirely dependent on just one function Gsik 2 ). Note that 
due to Lorentz invariance, is a function of the Lorentz invariant quantity namely, 
k 2 . When Lorentz invariance is broken in thermal field theory, one has rotational 
invariance only. For this situation, it is convenient to use the spatial projection 
operator P^ v which is defined as 

Pi = Pi = Pl = 0, = 6 V - ^. (3.8) 

But for 1 + 1 dimensions, it is easy to see that this tensor vanishes identically. 
Therefore, P^ defined in ( |3.4|) is purely longitudinal in this case. The only other 
projection operator in 1 + 1 dimensions which is symmetric and constructed out the 
momenta and r]^ is the 

*V = ^r- (3-9) 
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But tensors constructed from the above operator do not satisfy the Ward identity 
( |3.3|) . Thus when Lorentz invariance is broken, one only has the following change. 
The arguments of the functions in and Gs changes from the Lorentz invariant 
quantity k 2 to (a;, q). Thus the two point function can be written as 

Gwafiiw, q) = Pp,PapG B (w, q). (3.10) 

Writing it explicitly, we obtain 

Gu,tt = , 9 - — ^G B (u,q), G t t,tx — 7—5- — K^Gs(uJ,q), (3.11) 
(ur — q 1 ) 1 (uj z — q 1 ) 1 

Gtt,xx = 7—5 ^Gb{^i q), G tx ,tx = 7—5 ^Gb(uj, q), 

uj^q ui^ 

Gtx,xx — 7-5 ^Gb{&, q), G xxxx = - — 5 ov?G , q) ■ 

{u) z — q A y (uj z — q z y 

Thus all components of the thermal Green's function of the stress tensor are deter- 
mined by a single function Gb- 

3.2 Poles in the correlators 

We show here that the function Gb, which determines the thermal Green's function 
must exhibit a hydrodynamic singularity due to the propagation of sound. Using 
Lorentz invariance, the stress tensor of a fluid in 1 + 1 dimensions is given by Jl9| 

= (e + P)uV + Pif v - i{u»u v + r]^)d x u\ (3.12) 

where uF" is the 2- velocity with = — 1 and £ is the bulk viscosity. To obtain the 
linearized hydrodynamic equations, consider small fluctuations from the rest frame 
of the fluid. We then have the following 

T 00 = e + 5T 00 , T ox = 6T ox^ T xx = p + 6T xx^ ^ 13 ) 
u° = 1, u x = Su x . 

Note that u° = 1 up to the linear order due to the constraint = — 1. From the 
form of the stress tensor in ( p. 12 ), we can obtain the spatial variation of the velocity 
in terms of the stress tensor to linear order 

dj u * = (3.14) 
e + P V ; 

Substituting this for the velocity in T xx , we obtain to the linear order 

T xx = P + 5T XX , (3.15) 
= P ^dJT 0x . 



e + P 
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The hydro dynamic equations are <9 M T^ = 0. These reduce to 

d T 00 + d x T x0 = 0, d T 0x + d x T xx = 0. (3.16) 



Substituting the form of the linearized form of the stress tensor given in ( |3.15| ), we 
obtain 



d 5T w + dJT" x = 0, (3.17) 

doS T 0x + ^dJT 00 - -^-d 2 JT 0x = 0. 
oe e + P 

Here we have used the fact that the thermodynamic variable P depends only on the 
energy density e by some equation of state. Performing the Fourier transform of 
the above equations in position and time, we obtain the following algebraic set of 
equations for 5T 00 and 5T 0x 



oo5T m + q5T" x = 0, (3.18) 



-icuST 0x + iqv 2 s 5T 00 + q 2 -!L-5T 0x = 0, 

e + P 

where we have defined the speed of sound v s as 

o OP . 
vl = — . (3.19) 

Eliminating 5T 00 using the first equation, we obtain the following equation for fluc- 
tuations in ST 0x . 

,2 , ■ 2„,2 i S , 2 \ c- r 0.r 



-luj 1 + tq 2 vj + j^-pooq'j 5T 0x = 0. (3.20) 
Therefore the fluctuation in 5T 0x obey the dispersion with u given by 

uJ 2 -v 2 s q 2 + t-j-j 5 u J q 2 = 0. (3.21) 

Solving this to the leading order, we obtain the following dispersion relation for this 
longitudinal mode which we call the sound mode. 

u = ±v s q - i—-^ -q 2 . (3.22) 



From the equations in (|3.15|) and ( 3.18 ), it can be seen that the remaining fluctuations 



also obey the same dispersion relation. This implies by the usual arguments of linear 
response theory [[2(J that the two point function of the components of the stress 
tensor has a pole at the above value of u. Thus we find the function Gb in the 
retarded correlation functions has a pole at Q3.22 ). 
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4. The sound channel in gravity 



As discussed in the earlier section, the Dl-brane theory admits a holographic gravity 
dual for the ranges in temperatures given in ( p,12| ) and (p. 14 ). In this section, 



we first study the gravity solution of the Dl-brane and isolate a diffeomorphism 
invariant perturbation which we identify as the longitudinal mode corresponding 
to the sound mode. From the analysis of its equation of motion and by imposing 
quasi-normal mode boundary conditions, we derive its dispersion relation. We then 
evaluate the two point function of the components of the stress tensor of the Dl- 
brane theory holographically and confirm that it has the structure predicted by the 
general properties of hydrodynamics in previous section. 

To simplify our analysis, we first consistently truncate the 10 dimensional near 
horizon geometry of the Dl-brane in ( |2.8|) to 3 dimensions by dimensionally reducing 
on the 7-sphere using the following ansatz. 

ds 2 10 = e- UBir) g^(x)dx»dx u + e 2B{r) L 2 dn 2 , (4.1) 
= e -UB(r) (_4( r ) d f2 + c \{r)dz 2 + c\dr 2 ) + e 2B ^L 2 dS 2 x? . 

Using this ansatz in the 10-dimensional supergravity equations of motion, one obtains 
a set of coupled differential equations for the fields cr(r), Cx(r), cr{t), <p(r) and B{r). 
It can be shown that on identifying 

B(r) = ~4>{r) (4.2) 

and keeping the Ramond-Ramond flux through the 7-sphere constant, one can obtain 
a consistent truncation of the 10-dimensional equations to effectively 3-dimensions 



[jig , gj]] . The truncated set of equations of motion can be obtained from the following 
Einstein-dilaton system in 3 dimensions with action 

P. 



R-'-d^d^-V^) 



(4.3) 



Here, (3 — ^ and G3 is a three dimensional Newton's constant. The dilaton is denoted 
by and its potential is V = — -pes^. The coefficient in the dilaton potential is 
determined from the contributions due to the background flux through the 7-sphere 
and from its curvature. The equations of motion are 

□* - £M. (4.4, 

The Dl-brane in 10-dimension given in ( |2.8| ) reduces to 

ds 2 = -c T (r) 2 dt 2 + c x (r) 2 dz 2 + c R (r) 2 dr 2 , (4.5) 
0=-31og(£), (4.6) 
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with the components of the metric given by 



Cj-< 



with / = 1 — For future reference, we write down the the equations of motion 
explicitly in terms of these functions: 

c x cx cr 4 2 

^-^ + ^ /2 + 4? = 0, (4.9) 
4p + 2 ^-^/ 2 = 0, (4.10) 



c T c x 2 



(4.11) 



V. cr J P 

4.1 Linearized equations of motion for the perturbations 

We consider very small wave like perturbations in the background of the above so- 
lution g^y — ► g^ v + 5g fMU and 0-^0 + 5(f>. Due to translational invariance along the 
Dl-brane directions, we can assume that all perturbations can be expanded using its 
Fourier mode. Focusing on one such mode, we have 

6g^(t,z,r) = e-^-^V(r), 8<f>{t,z,r) = e^^^r). (4.12) 

We further parametrize the metric perturbations as 

hu = c T H a , htz = c x H tz , h zz = c 2 x H zz . (4.13) 

Following 0, [| we fix the gauge by choosing 5g rfl = 0. The equations of motion 
up to linear order in perturbations are 



Hi'+ln' 



\ctCrJ 

H 'L + In' (^) H' zz - H' tt ln'(c,) + °\z, + 2c||^ 

<p" + In' + CR (4-ir)<P+ - K) 

Cr \c t c x J l 



c\d^V_ 

(3 d(f) 2 



0. 


(4.14) 


o. 


(4.15) 


o, 


(4.16) 


o. 


(4.17) 



where A t = q 2 ^r and Z Q = A t H tt + 2quH tz + u 2 H zz . We also obtain following first 

c x 

order constraints from Einstein equations for 5R rfl . 

H' zz + 2-H' u - — 2 A t In' ^H tt + — 2 In' ^Z + (3^ = 0, (4.18) 
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H' tt - In' H tt + ^H' tz - M V = 0, 

\c T J qc T 



(4.19) 



r 2 &P 

W(c T )H' zz - \n\c x )H' tt - /30V + "f^o + c 2 R — V = 0. (4.20) 



In the appendix A, we show that the above 3 constaints can be consistently imposed 
on the 4 dynamical equations of motion. 

4.2 Diffeomorphism invariant sound mode 

Fixing the gauge 5g^ r = does not exhaust all gauge degrees of freedom. One is left 
with the residual gauge freedom under the infinitesimal diffeomorphisms x^ 1 — *> x^+^ 
with fi G {t, z, r}. Following the approach of [22|, we will construct a quantity 
using the above perturbations which will be invariant under the diffeomorphism 
of the metric i.e. 5g^ u — > 5g^ — V M £ V — V^. The perturbations change under 
diffeomorphism as 



c 



t)' 



^ - - 2ig^ - ^|^ r . (4.21) 

Then, the combination Zq = A t H u + 2quH tz + w 2 H zz changes as 

q 2 q 2 (4)' 2gcu,. w 2 / . f (4)' 



Z -> Z + — (2iu;)& + -g 2~fr + -^-(^6 - ^6) - — I 2zg^ H — £ 

C X C X C R C X C X V C R 

5Z = 2 ME^M (4.22) 



where 



A H = A^-u 2 . (4.23) 
In c x 



The dilaton also changes under diffeomorphism as 



- - ^£ r . (4.24) 



We find the following combination gauge invariant. 

Zp = Z + A v ip where A y = 2Ah ^}° x) . (4.25) 

Note that unlike the case of higher dimensional Dp-branes withp > 2, studied by [JIS], 
there is only a single gauge invariant mode for p — 1. Note that the Z constructed 
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Ill 



15 1 for p > 2 cannot be trivially extended for this case as H as defined in [15| 



does not exists for p = 1. Here, we note that the role of H for p — 1 is played by the 
dilaton fluctuation. We will call the diffeomorphism invariant fluctuation Z p as the 
sound mode as it constitutes fluctuations longitudinal to the wave directions. 

Next we outline the steps involved in obtaining the second order equation satis- 
fied by this mode. From the definition, we have 



Z" P = A"H U + 2A' t H' tt + A'^cp + 2A>^' 

+AH" t + 2quH? z + u?H" zz + Aptp" 



where 



A, 



2 c T 



-x 



(4.26) 



(4.27) 



Using equations ( ^.14[ )-( f4~2Ti| ), we simplify above to 



f(A t -LU 2 )Z P 



-4A t 



ln'(^ 

cx 



4AlnM — 

cx 



if + lh 



+ A" + 2c\V\A t - u 2 ) + -±V"A 







+2A'^'. 
Since A' t = 2A t W (g), 



(4.28) 



A", = 4A 



In' 



or 



1 2 



+ 2 A t In" 



Ct 

c x 



Evaluating the difference of the Einstein equations 
Using equation ( |4.29| ) and (|4.30|) , we obtain 



we obtain 



ct 
cx 



(4.29) 



(4.30) 



Z P + ln' ( C -^-)z' P - C 4(A t -u 2 )Z P 

V C R J C T 



i.bi..': • 



+2A'<p 



ct 
c x 



<J + In' ( 



CtCx\ 



V C R J 



A'+A"+2c 2 R V\A t 



u z ) + ^V"A V 

P 



Using Einstein equations and the dilaton equation of motion ( |4.8| )-( |4~Tl 
proove the following identities (see Appendix. O ). 



(4.31) 
, one can 



^\A t - lu 2 ) + A' = - 



+ ln' 



(-) 

\ctCrJ 



A v , (4.32) 



13 



AA t In' ( ^ 
1 c x 



</ + In' f 



\ C R J 



A' +A" + 2c\V\A t 



UJ 



2A' 



(4.33) 



Using relation ([4.32 ), one can write the relation ( |A.4| ) as 



A u 



CR 



t)<p + <p' 



Z;, \u ! [^)Z r 



(4.34) 



Finally using relation (|4.33|) , the equation (|4.31|) can be written as 



Z » + In' (Wx) Z' P - 4(A - w 2 )Z P = 2A' V 

\ CR J Cj. 



(4.35) 



Using the relation ( [4.34 ), we obtain 



Z'' + 



In' - 2^ 



R 



{A t -u 2 ) + 2 



A' 



A, 



= 0. (4.36) 



This is the equation for Zp from which we will obtain the dispersion relation for its 
quasi-normal mode. Here we make the following observation regarding the equation 
of Zp. Consider a minimally coupled scalar \l/ in the background ( |4.5|) , its equation 
of motion is given by 



1 



-A W^ggTdv*) = 0, 



(4.37) 



< + ln'(^*V-4 



(A t -u 2 )4> = 0, 



where g^ v is the background metric in ( |4.5| ). The second line is obtained from the 
first by focusing on the Fourier component \l/(t, z,r) = e ~ l ( UJt + c i z ) . Now comparing 
the above equation with the equation for Zp we see that the the equation for Zp 
( [4.361) is almost a equation for a minimally coupled scalar except for the terms that 
are proportional to A'/A^. Note that the dilaton dependence of the equation for Z P 
is entirely contained in these terms. 

Finally we emphasize that the analysis to obtain the equation for Zp in ( [4.36] ) 
depended only on the fact that the background is of the form (|4.5|) with a radially 
dependent dilaton profile (|4.6[) that satisfy the equations of motion ( |4.8D , ([4.9|) , Q4.10| ) 
and (|4.11|) obtained from the Lagrangian in (|4.3j ). At each step we have not assumed 
any form for the functions ct, Cx, cr, 4>, the dilaton potential V and any specific value 
for the constant f3. Therefore for any radially symmetric solution of the form (|4.5|) 
and a dilaton profile which solves the dilaton equation of motion, the equation for 
Zp is given by (|4.36j ). 
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4.3 The dispersion relation of the sound mode 



For further analysis, we will use the explicit expressions of metric coefficients and 
dilaton. We first change the variable to Y = and the independent variable to 



u 



j£, then the equation becomes 



u 



u{l — u 3 ) 



d u Y 



+ 



18m 4 (4A - 3) 



q 2 L 6 (1 - A - u 3 

Ir? (1 -u 3 ) 2 ' (l-u 3 )(4-4A-u 3 ) 2 



Y = 0, (4.38) 



where A = First we look for its solution near the horizon. For 1 — u <C 1, the 
equation becomes 



$y- 



l-u 



d u Y + 



a 2 oj 2 



+ 



6 



9(1 -u) 2 (3-4A)(l-u). 



y = 0, 



(4.39) 



where a 



L 3 
27^ 



. To pick up the ingoing solution at the horizon, let us define x 



ln(l — u). Then in terms of x and for 1 — u « 1, the above equation reduces to the 
oscillator equation 

F = 0. (4.40) 



d 2 + ^ 
x ^ 9 



Therefore, the ingoing solution for 1 — u « 1 behaves as e s c 
next consider a solution of the type 

Y = {l-u z )-^Z(u). 

Then the equation for Z(u) is given by 
d u Z 



We 



dtz 



u(l — u 3 
+Z 



-[2 + u 6 - 2iauju 3 



2 2 

a uj 



u 



a 2 q 2 



18w 4 (4A - 3) 



(1-m 3 ) 2 (1-m 3 ) (1-m 3 )(4-4A-m 3 ) 2 



0. (4.41) 



In this equation, all terms are dimensionless. From the definition of temperature in 
(|2.10|), we see that a ~ 1/T, therefore in the hydrodynamic limit 



uj « T and q « T, 



(4.42) 



we ignore terms of order q 2 /T 2 , u 2 /T 2 , uq/T 2 and higher, but keep terms of order 
lu/T, q/T . Performing this limit in the equation for Z, we obtain 



BIZ - {2 + (l-W} a , z 



18w 4 (4A - 3) 



0. 



(4.43) 



u(l-u 3 ) (1-m 3 )(4-4A-m 3 ) 2 
The well behaved solution of the above equation at u — 1, the horizon, is given by 

6A - 2(1 - A)(3 - Uauo) - w 3 (3 + 2iau) 



12(3-2W)(4-4A + u 3 ) 



(4.44) 
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To obtain the quasi-normal modes of this solution, we need to impose the Dirichlet 
condition Z = at the boundary at u — > oo. This leads to the following cubic 
equation for u. 



-4iau 3 + 6u 2 + Aiaujq 2 



3g 2 



Solving this equation perturbatively by assuming uj 
dispersion relation for the sound mode. 



1 



a 



to 



±7 =g- V + ... 



= 0. (4.45) 
q, we obtain the following 

(4.46) 



with 



a 



z ' o 



4vrT 



4.4 Holographic evaluation of the stress tensor correlators 

In this section, we use the standard prescription of the gauge/gravity correspondence 
to evaluate the stress tensor correlations. For this, we first need to expand the 
action in (|4.3|) along with the Gibbons-Hawking boundary term to second order in 



(4.47) 



the fluctuation Hu V - The bulk action and the boundary term is given by 



S 
S 

where 



Sbulk + S, 
1 



GH, 



16ttG, 



d 3 Xy/—g 



2 ^ 



-v(<f>) 



d 2 xV —hK\, 



16 n 

— V = - 
9 L 2 

and h is the boundary metric at a large but fixed value of r and K is the extrinsic 
curvature of the boundary 4 . Using the equations of motion ( f4.14| )-( f4.17| ) and the 
constraints ( f4.18p - Q4.19p we can rewrite the expansion of the bulk action to second 
order in the fluctuations as a total derivative in the radial coordinate. This is given 
by 



S, 



(2) 
bulk 



+ 



d 3 x 
ctCx 

CtCx 

(3 c T c x 
4 c R 



34 



2c T c R 



H tzH' tz + 



-x 



CtCr \c t 



— - 2— | H, 



cx 



Hz 



(3 c T c x , 

(pip 

2 c R 



[H tt + H zz )(H tt + H zz )' + ^ W(c T c x )H tt H : 



2c 



R 



HttH' tt + H ZZ H' ZZ + (H tt — H z 
- i 

(H tt - H zz )(j)(p . 



c T 



z2L 
cx 



H z 



(4.48) 



4 In general there are counter terms one has to add to regulate the action (H, ^I], ^[ > these 
counter terms are crucial to regulate the one point function of the stress tensor. However for the 
two point functions which we will be interested in, they are not relevant. 
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The second order perturbation in extrinsic curvature term is given by 



-hK 



r 3 



8c T c R 



i ( °±Hi - 3 C ^Hl 
C T c x 



^H' tz H tz 



+ ^(H tt + H zz ) (2(H a + H zz )' + \n'(c T c x )(H tt + H zz )) 



(4.49) 



Since the bulk action at second order in perturbation is just a total derivative we 
just have to evaluate its contribution at the boundary. Then the boundary action at 
second order in fluctuations including the Gibbons Hawking term then reduces to 



1 



16ttG ? 



dudq 
12^ 



(4.50) 



where 



C 



-x 



-Zp{k)H' tz + /^contact) 



Auqc T c R 

and /Contact represents the part of the Lagrangian without any derivatives. 

H tz H t t 



(4.51) 



C T C X 



Ac 



R 



ln'(c T c x )H tt H zz + °-^Hl + C ^H 2 ZZ - 2^ In' 
cx ct uj 



cx 



C X C> Xtt2 CtCx 
H n l? : if 



<-]■<■ H <X ' 4c R 
c_x_ 

c T 



2ip 



c\V 



pin'(c T c x )} -2p<t>\H tt - H zz ) 



H tt + 



2Z \c\ In or ( c x 



In' 



Ct 



(4.52) 



To obtain the above form of the action we have used the constraints ( |4.18|) -( f4.20|) 
to rewrite all derivative terms in C in terms of the derivative H' tz . We can further 
manipulate the boundary action and reduce it to the form 

1 



16ttG, 



%%A{u>, q, r)Z' P (r, k)Z P (r, -k) + 5, (2) 
(2ny 



CTi 



where 



A(u,q,r] 

The contact term in equation ( |4.53 ) is 

dujdq 



P ctc x 
'2AI c R • 



(4.53) 
(4.54) 



o(2) _ 1 



(27T) 



A 



A' t H tt 



2A t + A * 



|4 + MiA* - uj 2 ) - A v ( m'(c K ) - r_ 
+Z (-ln' 



c_x_ 

ct 



A r 2 A 

->4~m^ (ct) 



cx_ 

ct 



In' 



Cx 



ct 



(4.55) 
Z P (r,-k). 



Now that we have the boundary action, we can evaluate the stress-tensor cor- 
relators by using the standard rules of the AdS/CFT correspondence. Note that in 
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this case the geometry is not asymptotically AdS. However we expect the rules of 
the AdS/CFT correspondence to still be valid for this case since it can be related to 
the anti-de Sitter geometry up to a conformal factor ]T^, [14]]. The boundary values 
of the fluctuation H^ u couple to the stress tensor of the boundary theory as in || 

Coupling = ~ J d'x\HlT* + H° ZZ T ZZ + 2#°T'*]. (4.56) 

Here, the indices are raised, lowered and contracted using the flat metric ds 2 = 
—dt 2 + dz 2 . The superscript indicates the fact that we are looking at the r — > oo 

r 2 

or the u = — ► limit of the corresponding bulk fields. Thus these correspond to 
the boundary values of the fluctuations H^ u . Using the above coupling and the rules 
of the AdS/CFT correspondence summarized in the equation 

(exp(zS coupling )) = e W [iSW(Hp], (4.57) 

we can evaluate the various two point functions of the components of the stress 
tensor. Consider the retarded two point function G tt ,tt = —i([Tu,T t t\), using the 
AdS/CFT prescription we obtain 

G tttt = -A =; (4.58) 

6Hg(k)6H&-k) 

From the definition of Zp, we note that its boundary value is related to the boundary 
values of the fluctuations by 

Z% = q 2 H° t + 2ujqH? z + uj 2 H° zz + A°y. (4.59) 

From Q4.44D , we see that the expansion of Zp near the boundary is given by 



Z P = C(l + ...) + Du 3 {l + ...) 



Z P 



Dr, 



1 + + 1 + 



Cr 



6 



(4.60) 



Here C(u, q) and D(u, q) are independent of u. The ellipses denote higher powers in 

2 

u = hj. Now substituting the above expansion for Zp, in it is easy to see that 
relevant term in the boundary action is the first term in ([4.53 ) which involves the 
derivative of Zp. This results in the following expression for the two point function 
of the stress tensor. 

G - ! ^ (2) 1 I' K 6 ^ (A61) 

tt ' U ~ SH&fySH&i-k) ~ 16ttG 3 (uj 2 - q 2 ) 2 U C 1 ' 1 

Comparing the above holographic result with ( |3.11| ) we see that the expected Lorentz 
factor for this correlator is reproduced. Furthermore we can read out the holographic 
value of Gp as 

1 6r 6 D 

G ^> = -issuer (4 - 62) 
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The poles in the Green's function are therefore same as the zeros of the factor C. 
The Dirichlet boundary condition for the mode Zp = at the horizon is equivalent 
to setting C = as noted for the case of backgrounds asymptotic to AdS which were 
studied in [f22|. Thus the poles in the two point function of the stress tensor which 



leads to the dispersion relation of the sound mode is as given in equation ( |4.46| ) . As a 
consistency check of our calculations, we have evaluated all the remaining two point 
functions of the components of the stress tensor. In each case the expected Lorentz 
factor given in fl3.11| ) is reproduced with the same expression for Gb{oj, q). 



5. The ratio of bulk viscosity to entropy density 

In this section, we use two methods to evaluate the ratio of bulk viscosity to entropy 
density £/s of the Dl-brane theory in the gravity regime which is valid for temper- 
atures in the range VXN-I « T « V\. We use the dispersion relation of the 
sound mode to read out this ratio, then as a cross check we determine this ratio by 
evaluating the bulk viscosity directly using the Kubo's formula for the bulk viscosity 
in one spatial dimensions. We then show that the ratio £/s continues to be the same 
for the temperatures in the range V'AiV" 1 << T << \^XN~s when the Dl-brane 
background is replaced by the Fl-string solution. In both these regimes, the ratio 

i = j_ 

s Arc 

We also show that for gauge theories corresponding to Dl-branes at cones over Sasaki- 
Einstein 7-manifolds the ratio continues to be 1/A%. 

5.1 £/s using dispersion relation of the sound mode 

Using the dispersion relation of the sound mode in gravity ( |4.46| ) and comparing 



it with the dispersion relation of the sound mode using general hydrodynamic con- 
sideration ( |3.21[ ), we deduce the transport properties of the Dl-brane matter at 
temperature T. As we have seen in section 3, general hydrodynamics considerations 
give the following dispersion relation for the sound mode 

% 1 2 2 dP 

Here, v s is the speed of sound in the medium P, the pressure and e the energy density. 
Now comparing this with the dispersion relation obtained from gravity in ( 4.46| ), we 
can read out the speed of sound in the Dl-brane matter as 

vl = \ (5.2) 

It is now easy to see that from the definition of the speed of sound v s = we 
obtain the equation of state 

e = 2P 
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for the Dl-brane matter. Thus the mediums seem to behave as a conformal fluid in 2 
spatial dimensions just like that of the M2-brane. The fact that the thermodynamic 
properties of the Dl-brane theory is similar to the M2-brane was noted earlier in |^6| 



and [37]. In [26], it was noted that the entropy density of the Dl-brane behaves like 
that of the M2-brane 5 , while in [37| the speed of sound in the Dl-brane matter was 
evaluated holographically just using thermodynamics of the Dl-brane and shown to 
be the same as given in Q5.2| ). It will be interesting to understand this coincidence 
since the Dl-brane gravity solution is certainly not of the M2-brane form. 

From comparing the dissipative part of the dispersion relation in ( 4.46Q and (|5.1| ) 
we see that 

^ = - = -— (53) 
e + P 3 6r 2 4vrT' 1 ! 

In the absence of chemical potentials, we have the following thermodynamic relation 



e + P = Ts. 

Substituting the above equation in ( |5.3|) , we find that 

i = L 

s An 



(5.4) 



(5.5) 



Using ( |2.10|) and the definition of L from (|2.2| ), the entropy density can be written 
in terms of the field theory variables and temperature as 



2 4 vr§ N 2 T 2 
3 3 V\ ' 



Then the bulk viscosity is given by 



(5.6) 



2 6 7ri N 2 T 2 
3 3 v^A ' 



(5.7) 



Let us now compare the result in ( |5.5| ) with the ratio of bulk viscosity to entropy 



density for Dp-branes with p > 2. Using the results of [T5|] for Dp-branes with p > 2 
we have 

£ _ £ri _ 1 2(3-p) 2 
s rj s An p(9 — p) 



(5.J 



Note that the above equation is valid p > 2, since the method used by [15j is valid 
only for p > 2. But using the result (|5.5| ), it is clear that the formula for £/s obtained 
for p > 2 continues to hold for also p = 1. 



3 See below equation 2.17 in 126 
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5.2 £/s using the Kubo's formula 

As a cross check of our calculations, we use the Kubo's formula to obtain the ratio 
of bulk viscosity to entropy density for the Dl-branes. In 1 + 1 dimensions, Kubo's 
formula for bulk viscosity is given by (see for instance in ||20|| ) 



£ = lim - ^ dt [ dze iu,t ([T zz (x),T zz (0)}). (5.9) 
= lim -G zz ^ zz (uj,q = 0). 

As we have seen in section 4.4, the two point function of the stress tensor from gravity 
is given by 

uj a 1 r 6 D 

G zz , zz = {uj2 _ q2)2 G B (u, q), G B = -^-JL- (5-10) 

Here, the coefficients C and D are related to the asymptotic expansion of the solution 
near the boundary at u —>■ given by 

Z P = C(1 + ...) + Du 3 (l + ...). (5.11) 

For the sound mode in fl4.44|) , we find that 

{-Mau 3 + Aiauq 2 + Qu 2 - 3q 2 ) 



C 



9(2iacu - 3) 

ufiiuq 2 + %ia 2 uu{u 2 - q 2 ) 2 - 12a(2g 4 - 3g 2 cu 2 + w 4 )] 



D = — 1 -, w ON 1 -. (5.12) 

54(3i + 2acj)(g 2 - uo 2 ) v ; 

In limit q — > 0, we obtain 

c^Jf (, 13) 

and hence 

G zz , zz (u,0) = G B (u,0) = (5.14) 

Now using the temperature and the entropy density of the Dl-branes solution given 
in ( |2.1(J| ), we can write the two point function as 

G zz>zz (u,0) = -Ts^ = -i^. (5.15) 

Substituting this result for the two point function into Kubo's formula in (|5.9| ), we 
obtain the ratio 

- = 7~- ( 5 - 16 ) 

This agrees with that obtained from examining the dissipative part of the sound 
mode in (|5.5|) and therefore is a consistency check on our calculations. 
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5.3 Universality of the ratio £/s 



(i) £/s for the fundamental string 



We have seen that in the range of temperatures y/XN 3 << T << \f\ for 
which the Dl-brane solution can be trusted, the ratio £/s is given by the ideal value 
1/4-7T. We now show that the ratio continues to be 1/Att for the temperature range 
« T « 3 for which the fundamental string solution can be trusted. 

From ( |2.13| ), we see that the only difference between the Fl-string solution and that 
of the Dl-brane solution is that of the dilaton profile. In fact the Fl-string solution 
can be obtained from the Dl-brane solution by replacing <fi — > —cf). On compactifying 
to 3 dimensions just as for the Dl-brane case, one obtains a consistent truncation of 
the equations of motion. The equations of motion can be derived from the action in 
3 dimensions given in (|4.3|) but with 



16 



24 4 
— p 

L 2 



(3 = — , and V 
The Fl-string solution in 10 dimension given in fl2.13|) reduces to 



ds 2 = -c T (r) 2 dt 2 + c x (r) 2 dz 2 + c R (r) 2 dr 2 
■ r ■ 



3 log 



L 



(5.17) 



(5.18) 
(5.19) 



with ct, cx, cr and / given by Q4.7Q just as for the Dl-brane case. Note that the only 
change from the Dl-brane case is that of the sign of the dilaton. 

The analysis of the linearized perturbation and obtaining the equation for the 
sound mode in section 4.1 and 4.2 just depended on the fact that the background is 
of the form in ( |5.18 ) and radially symmetric. Thus the equation for the sound mode 
for the Fl-string solution continues to be ( |5.20| ) which is given by 



Z P + 



Kr) 


A' ' 
- 2—^ 


z' P - 




_ 





r 2 A' 
^(^-^ 2 ) + 2-^ln' 



(■-> 



A 



■p 



cx 
c T 



0. (5.20) 



In the above equation, the only dependence on the dilaton profile is due to the ratio 



-j^. From the definition of A^ in ( 4.25|) , we see that this ratio is given by 



A; _ [A H \n'(c x )]> 
A v A H \v!(cx) 



(5.21) 



Since Ah just depends on the components of the metric Q4.23j ), the only difference 
for the Fl-string can arise from the ratio <fi" /(£>'■ But it is easy to see from ( |5.18| ), 

A' 

that this ratio is same as that of the Dl-brane case. Therefore the ratio -f- for the 
Fl-string remains identical to that of the Dl-brane case. Thus, the equation of the 



- 22 - 



sound mode Zp is the same for the Fl-string. Now the rest of the analysis to extract 
the dispersion relation for the sound mode proceeds identical to that of the Dl-brane 
case. This results in the ratio 

£ = — 

s Air 

for the Fl-string background. We have thus extended the evaluation of the bulk 
viscosity to the temperature range VXN^ 1 << T << \/\ for the U(N) gauge theory 
with 16 supercharges in 1 + 1 dimensions. 

(ii) £/s for Dl-branes at cones over Sasaki- Einstein 7-manifolds 

Here we show that the ratio of the bulk viscosity to entropy density for theories 
of Dl-branes at a cone over Sasaki-Einstein 7-manifolds is also given by l/An in the 
regimes given by ( |1.2|) . The near horizon geometry of the corresponding supergravity 
solution is given by 

ds 2 = F-t(r) (-f(r)dt 2 + dx{) + H$(r) (j^ + r 2 dS 2 x )j , (5.22) 
e 4>(r) = H{r)^, 
F 7 = 6L 6 ux 7 - 



where 



H(r)=(^\\ f(r) = l-(^y. (5.23) 



and dS x? stands for the metric of the 7 dimensional Sasaki-Einstein manifold X 7 , 
and u>x 7 is its volume form. Note that the only way this background differs from 
Dl-brane at flat space is by the replacement of the 7-sphere by the 7-dimensional 
Sasaki-Einstein manifold. To be explicit here, we write down the metrics of two 
7-dimension Sasaki-Einstein manifolds known as Q l ' 1)l and M 1,1 ' 1 in the literature. 
The metric of Q 1,1,1 is given by 



dSLi,! = - V cos^^) 2 + - V (d^ 2 + sin 2 0,^), (5.24) 

v 16 z — ' 8 z — ' 

i=l i=l 



while the metric of M 1,1,1 is given by p9| 



where 



dS'Li.i,! = — (dr + 3 sin 2 /icr 3 + 2 cos Ozdfa) 2 (5.25) 
64 

3 X o 1-2 / 2 2 2 2* 

+ -\djj + - sm fj,(cr 1 + cr 2 + cos /xcr 3/ 
+^ 2 + sin 2 9 2 d</%), 



(Ti = d9i, <j 2 = sin^id^i, a 3 = (dip + cos^d^). (5.26) 
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Similar to the case of flat space we now dimensionally reduce to 3 dimensions using 
the ansatz 



ds 2 l0 = e- UB ^g^(x)dx^dx v + e 2B(r) L 2 dS 2 x? , (5.27) 
Again one can verify just as in the case of flat space the identification 



B(r) = ~4»{r), (5.28) 

provides a consistent reduction provided the 7-form flux through the Sasaki-Einstein 
manifold is held constant. The equations of motion reduce to the Einstein-dilaton 
system in 3 dimensions. The coefficient of the potential of the dilaton in this system 
is dependent on the the Ricci-scalar of the Sasaki-Einstein manifold. For the metric 
L 2 dS Xr when X 7 is either Q 1 ' 1 ' 1 or M 1 ' 1,1 the Ricci-scalar is given by 

R(X 7 ) = H. (5.29) 

This is same as that of the 7-sphere, which implies the coefficient of the dilaton 
potential will remain the same in the effective action. Note that the Ricci-scalar for 
the general Y Sasaki-Einstein 7-manifolds constructed in |30[ and recently studied 
in 

in 3-dimensions is given by 

1 f „ , / , s 8„ , 24 



31] is also given by ( |5.29 ) 6 . Taking all this into consideration, the effective action 

d 3 x^9 (r( 9 ) - h^(j> + ^el^j , (5.30) 



16wG 
where 



1 L 7 Vol(X 



7) 



r c ^ (5 ' 31) 

We see apart from the definition of the 3 dimensional Newton's constant, the effective 
action remains the same as that of Dl-branes in flat space. The background values of 
cx{r), cx{r), Cr(t) and 0(r) is also same as that of Dl-branes in flat space. Therefore 
the quasi-normal mode analysis remains identical and the ratio of the bulk viscosity 
to the entropy density remains the same. The analysis for the case of the Fl-strings 
at cones over Sasaki-Einstein 7-manifolds is also identical to that of Fl-strings in flat 
space and yields the same ratio of bulk viscosity by entropy density. 

The gauge theory of Dl-branes at cones over Sasaki-Einstein 7-manifolds is dif- 
ferent from that of the SU (N) theory with 16 supercharges for the case of Dl-branes 
in flat space. The number of supersymmetries and the matter content of the theory 
is different. Inspite of this for the temperature range ( |1.2| ), the theories of Dl-branes 
at cones over Sasaki-Einstein 7-manifolds have the same value of £/s = 1/47T as that 
of Dl-branes in flat space. 



3 See below equation 2.6 of |31[ 
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6. Conclusions 

We have shown that the ratio of bulk viscosity to the entropy density for the SU(N) 
gauge theory with 16 supercharges in 1 + 1 dimensions on the Dl-branes in the 
temperature range is given by 1/Att. For temperatures 

outside this range, the Dl-brane gauge theory flows to a conformal field theory. We 
therefore expect the bulk viscosity to entropy density of this theory to vanish for 
T « and T » vA. A technical result of our analysis is the equation for 

the sound mode given in (j4.36|) for any Einstein-dilaton system in 3 dimensions given 
by the Lagrangian ( f4.3|) which admits a radially symmetric solution of the form (|4.5| ) 
and a dilaton profile determined by the dilaton equation of motion. 

We have also seen that for the theory of Dl-branes at cones over Sasaki-Einstein 
7-manifolds the near horizon geometry dimensionally reduced to 3-dimensions is 
identical to that of Dl-branes in flat space. This implies that the ratio of bulk 
viscosity to entropy density is given by 1/Att. At this point, it is perhaps worthwhile 
to investigate other 1+1 dimensional systems which admit gravity duals with different 
near horizon geometry. There are two possible interesting geometries one could 
explore: One can turn on a R-charge along the Cartan directions of the SO (8) R- 
symmetry of the Dl-brane theory. This results in a R-charged black hole in the 
Einstein-dilaton-Maxwell system in 3-dimensions. The near horizon geometry of 
such black holes in 3-dimensions is different and it will be interesting to evaluate 
the ratio of bulk viscosity to entropy density in this geometry. Another interesting 
geometry in 3-dimensions is that of the D1-D5 system. This geometry is conformal, 
however one could turn on a relevant operator in the orbifold theory of the D1-D5 
system which will render it non-conformal. The holographic dual geometry to such 
a system will necessarily have a non-trivial dilaton profile and therefore a non-trivial 
bulk viscosity. This system will fall into the general Einstein-dilaton system studied 
in this paper, with an specific dilaton potential. We expect the dual geometry to 
be radially symmetric since one can possibly choose the relevant deformation to be 
translationally invariant along the brane directions. From the analysis in this paper 
we can then conclude that the equation for the sound mode is given by ( }4.36|) . It 
is tempting to speculate that the ratio of bulk viscosity to entropy density for this 
system will also by 1/Att. It will be interesting to construct this dual geometry for 
the deformed D1-D5 system explicitly and verify this conjecture. 

One of our motivations to study holographic duals to hydrodynamics in 1 + 1 
dimensions is the possibility of taking the non-relativisitc limit to obtain the Burger's 
equation 7 and study turbulence in one spatial dimension. However the scaling put 
forward in |32|, |33|] to obtain the non-relativisitc limit results in velocity of sound 

7 Burger's equation is the non-relativistic Navier-Stokes equation in one spatial dimension given 
by d t v + vd x v = £,d%v, where v is the velocity. 
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being infinite. This leads to hydrodynamics with divergenceless velocity flow, 



drf = 0. 

In one spatial dimension, this implies a trivial solution for the velocity field. The 
velocity field is constant in both space and time. It will be interesting to explore 
the possibility of taking other non-relativistic limits which can lead to the Burger's 
equation. 



Note added: After completion of this work, we noticed the preprint |34| which has 
some overlap with this paper. 
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A. Consistency of the constraint equations 



Note that the equations ( }4.18| )-( f4.20| ) form 3 constraints for the 4 dynamical equa- 
tions ( [4.14[ )-( P:.17| ). These constraints must therefore be consistent with the dynam- 
ical equations. That is, on evolving the constraints by the dynamical equations, one 
should not generate new constraints. In this appendix, we show that by differen- 
tiating the 3 constraints in ( 4.18 ) to ( 4.20|) and reducing them again to first order 
equations using the dynamical equations in ( |4.14j ) to ( |4.17|) , one does not generate 
any new constraints. This allows us to conclude that the constraints are consistent 
with the dynamical equations and they reduce the number of dynamical variables to 
just one, which we have identified as the sound mode Zp. 

We start with writing down 2 equations which we will use in our analysis. We 
define 

Z = A t H tt + 2quH tz + u 2 H zz , (A.l) 



where A t = q ^f. Its first derivative is 



Z' = A' t H tt + A t H' tt + 2uqH> tz + to 2 H' zz . 



(A.2) 
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Using the constraint equations, one can remove H' tt and H' zz in favour of H' tz . Then 
the expression reduces to 

Z' a = j^ + ^ln' (^\ Htt -uj 2 \n' fe) H zz -2uq\n' fe) H tz + {A t -u 2 )(36'^. 

(A.3) 

Next we write down an equation obtained by differentiating Zp. 
Z' P = Z' + A^' + A'^, 

= - In' fe) Z P + jln' fe) A„ + (A t - u 2 )(36' + A'^ <p + A^'. (A A) 

Here we have used the equation ( |A.3| ). Equation ( A.4 ) can also be obtained by 
adding A t times ( |4.18[ ) and u 2 times ( |4.19| ). 

To check the consistency of the constraint equations, we first differentiate them, 
then substitute the dynamical equations and reduce the equation to an equation 
with at the most first derivatives. We then show that we do not obtain any new 
constraints. We start with the equation ( |4.19|) . We get 



-In' 



^) H' a + Kz In'(or) + ^ In' (-^-) K - In" ( C -*) H tt + %Z P 

Cr J qc T 2 \c x CtJ \c T J 4 

" J <p - (36'p' = 0. (A.5) 



O 2 &P C R 
+ I 2c R -^- 2" A(p — 

Crp J 

When we use ( |4.30| ) on the combination of equations ( |A.5| )-( f4T2"U| )+hi' f^r) times 
( gT|) , we get 

ttp = 0, (A.6) 

where 

, / c T c x \ 



tt = c 2 R V - In' 



(A.7) 



V cr. J 

Now ( |A.6| ) is consistent for any fluctuation p because tt = 0, since the back ground 
statisfies the dilaton equation of motion ( 4.11]) . Thus the equation ( |A.5 ) is consistent 
with equations of motion. 

Next we differentiate ( 4.181) and use the dynamical equations to remove the 
second derivatives. Doing this, we get the equation 



-In' 



c tc 2 x 
cr 



Hi 



Q 



-x 



In' 

LO \ C T C R 



)^ + {ln'(cx)-^ln' 



c T 



+ 



-^<U„lnM^ 



Or 



+ 



C X 



c T 



r 



k" C R a i dV 



cx_ 

c T 
cx 



HL 



c T 



+ { -— A v hi 



cx_ 

ct 



+ (36'\p' = 0. 



(A. 
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Note that using equation (|4.30 ), one can show that 



In" [ — 1 - 

c T 



ln'(^ 

c T 



l2 = ln' (°A ln'(£ 



cx 



CR 



(A.9) 



Using ( |A~9l ) on the combination of equations (WM - (]Oo|) +^ln' (fH times (gjj) 



+ln' ^ 



times ( |4.18| ), one obtains 



-Clcp = 0. 

Hence, the equation ( |A.8|) can be written as a linear combination of earlier 
equations. 

The derivative of third constraint equation (}4.20|) gives us 
ln"(cT) — In' Ct In' ( T x ) — In' Cx In' + ^-(f)' 2 

\ Cr I '1 



(A.10) 
constraint 



+H' tt 



- ln"(c x ) + In' c x In' f ^* + In' c x In' c T - § 
/In' ^"j-^-^ + c^' 1 



+4v 



~ ^2 



\ Cr J C T 

# [-2P- , n < (™) + ^ - 4) + ^ In' 4' 

L WW V4 4/ 4 V 4 / 4. 

|{^-ln'(^Z P }=0. (A.11) 

We use equation Q = to remove the /3(f)"<p' term and equation ( |A.4|) to 
Next we add 2 In' times flOOl) to it. We note that equation (p 

be written as 



■ remove Z' P . 



In" ct + In' ct In' ( T X ) = In" cx + In' cx In' ( T X ) . 

V C R J \ c R J 

Using this equation with the above manipulations, we obtain 



(A.12) 



In" c T + In' c T In 



- 2 In c T In c x + -<p 

V c R J 2 _ 



(H zz -H tt )' = 0. (A. 13) 



The term in the square bracket is a combination of Einstein equations of motion, 
namely equation (|4.9|) -^( fLlO| ), hence it vanishes. This implies that (|A.ll ) can be 
written as a linear combination of the three basic constraints upto equations of 
motion. 
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B. Identities from background field equations 



(i) Derivation of equation ( 4-32 ) 



First we will derive equation ( 4.32 ) from the background equations of motion. 
The (a; 2 ) and (q 2 ) part of the equation should independently vanish. Since 



A' — ——A + 1 



A^Wrr-l 21nV / ln'( ^ ) \-^\n"< x 



we need to show 



k'2 



2 
and 



In' 


m 






2 






2 



+ In" cx + In' Cx In' 



— 

ctCrJ 
cx \ 



CtCrJ 



Using equations (|4.9|) and (|4.8| ), we obtain 



ln"c T = ln'c T ln' 

\ctJ 4 2 
\n" c x = ln'cxln ( — J - -i 



/2 _ 2rj>. 



(B.l) 

(B.2) 
(B.3) 

(B.4) 
(B.5) 



Replacing In Cyfln cy) in equation |B.2| |([ |B~3] ]), we find that the left hand side 



vanishes due to equation (14. lOf) . 

(ii) Derivation of 



Next, we will derive relation ( |4.33| ). The (uj 2 ) and (q 2 ) part of the equation 
should be independently satisfied. They are 



c\V" 



G' + 2c R V + 2-|— In c x + G 



m > ln ' ( ^ ] + K ' + 2c |p' + H^! l n ' CT + # l n ' 



0, (B.6) 



4 \ 20" , N 

-^V +— K = 0. B.7 



where 



G 



K 



ln c x - — In c x 



ln" c T + 2 In' c T In' ( — ) - — ln'c r 



To show (|B.6|) , we first evaluate G'. 

G' = 



20" 2 



ln^cx-^ln'cx 



(B.8) 
(B.9) 

(B.10) 
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Using equations (4.8) and ( 4.11 ), we can write 



ln (3) c x - ^—r In' c x = In" c x In' ( — ^ + In' c x In" c T + ^77 In' cx In' f CtCx ^ 



2c R V>c> x c' R c R V ". Wcx _ c 2 rVWcr 







(B.11) 



Then the left hand side of equation 



LHS = Gin 



+ 



^' 2 , / Z' CTCx 
— m 



becomes 

+ In" c x In' 



+ ^ln'c,ln' 



k'2 



In 



V cr ) 



%%P' c'r c'x 

13(f)' c R c x 



+ ln"cxln' ( % ] + In' c x In" c T 



CR 



+ In' c x In" c T 



-cIvWcr 



R 



+— In c x ln c R - 



P<? CrC X 



(B.12) 



Using equations (|4.10 ) and fl4.ll ), one can show that 



(In" ct In' cx + In' cy In" cx) = — 



2-P P 



'2 



c\P In' c R + 



In' 



ctc x \ 

CR ) 



Then we obtain 
LHS 



2\n'c R 



-21n"cx + ^-ln'cx- 



^Yx~ 2CrV 



(B.13) 



(B.14) 



Using equation ( |B.5| ) to replace In" cx and equation ( [4.1 1|) to replace 4>", the Lif S 1 
then vanishes due to equation ( |4.10| ). 
To show (p.7|) , we first evaluate 



<h" 2 
K = -2—K+ — 



ln"c T <'ln'(^]+21n'(^)l > + ^^ + 



T 



c T 



cx 



c' T c' R 



k'2 



<-j:<-r 2 

' or 



In' 



ctc x 
cr 



-c R V'4>' - c\V In' c R + 2 In' c T In" f^j + ^ In' c T In' 



cx 
/3 



(B.15) 



Using equations (|4.11|) and fl4.9|) , we can write 
C ^'ln'c T ln' Q? = 20"^^ 



CT Cr 



-1 + e 

ct 



+ 



c|P 



(B.16) 
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Using it, we get 



+ ■ 



In' CT ( I. In' )- In" c x - 2 In' c T In' ( ™L 



+ ln"c T ln' -clV In' (c T c x ] 



(B.17) 



The left hand side of the equation flB.71) then becomes 



LHS = - 



In" c T In' f ^ + 2 In' f ^ In' c T In' f + 2/30' 2 In' 
-c\V \n(c T c x ) - In' c T In" c x - 2(W c T f In' f ^) 



Using equation ( |B.13| ), we get 



LHS = - 



In' f?E?j |41n"c T + 21n'c T ln' 



J2 



Using equation 



LHS 



we get 



In' 



U2 



-2 4 P + 21nVw(-^-)} 



+ ln <^ -c R V + — - - 2(ln c T ) 



\ °R J 



In' 



4 In' c T In' c x + 2 In' c T In' 



+ ln' {2 In' c T In' c x -2(ln'c T ) 2 } 

= = RHS, 

where we used equation ( |4.10| ) to get middle step. 



ct 
ex. 



(B.18) 



(B.19) 



(B.20) 
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